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ABSTRACT 
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chamber  pressure  and  thrust  during  the  transient  periods  of  buildup  and 
decay.  The  results  are  differential  equations  and  solutions  of  closed  form. 
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LIST  or  SYMBOLS 


Cross -sectional  area,  surface  area 

Coefficient 

Constant 

Specific  heat  at  constant  pressure 

Constant 

Thrust 

Mass  of  propellant 

Combustion  index 

Pressure 

Ignition  pressure 

Gas  constant 

Burning  rate 

Temperature 

Time 

Volume 

Velocity 

Nozzle  angle  in  divergent  part 
Ratio  of  specific  heats 
Correction  factor 
Density 

Independent  variable 


LIST  OF  SYMBOLS  (Continued) 


Discharge  function 
Mass  flow  rate 

Refers  to  ambient  condition 
Location  at  chamber 
Location  at  exit  or  separation 
Refers  to  propellant 
Location  at  throat 
Initial  condition 


INTRODUCTION 


In  determining  the  total  impulse  and  performance  values  of  solid 
propellant  rocket  motors  usually  the  thrust  buildup  period,  which  is  the 
time  from  ignition  up  to  a  certain  rated  thrust,  and  the  ..rust  decay  period, 
which  is  the  time  from  an  average  effective  thrust  down  to  zero,  are  not 
included.  The  total  impulse  can  be  accurately  calculated  by  integration  of 
the  thrust  over  the  total  operating  time.  For  calculating  the  thrust  during 
these  transient  phases  of  buildup  and  decay,  the  chamber  pressure  as  a 
function  of  time  has  to  be  evaluated  first.  During  these  periods  quasi¬ 
steady  flow  through  the  nozzle  is  assumed.  General  equations  will  be 
derived  which  are  differential  equations  applicable  for  cases  when  several 
parameters  are  variables,  as  area  of  flame  front  changing  with  time,  com¬ 
bustion  index  changing  over  certain  pressure  ranges,  average  ratio  of 
specific  heats  changing  with  expansion  ratio,  etc.  Then,  assumptions  are 
made  where  those  parameters  are  constant  and  solutions  of  closed  form 
are  derived. 
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CHAMBER  PRESSURE  BUILDUP 


General  Derivations 

Consider  a  general  case  of  burning  fuel  in  a  solid  propellant  rocket 
motor  where  the  combustion  gases  escape  through  the  nozzle. 


By  applying  the  law  of  conservation  of  mass  in  the  chamber  the  following 
equation  can  be  used 


dp  V 

- ; -  =  0)r  -  W 

dt  i 


(1) 


where  o>£  is  the  rate  of  fuel  consumption  and  Li  is  the  discharge  rate  through 
the  nozzle. 

The  burning  rate  of  a  solid  propellant  increases  strongly  with  an 
increase  in  pressure  under  which  the  propellant  burns.  The  burning  rate 
can  be  expressed  as 


n  /->v 

r  =  apc  (2) 

where  a  and  n  are  coefficients  which  depend  on  the  propellant  composition, 
on  temperature  and  on  operating  pressure.  Introducing  the  surface  area 
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of  the  flame  front  Af  and  the  density  of  the  propellant  p^  the  rate  of  fuel 
consumption  can  be  written 


iif  =  A£pfapc 


n 


(3) 


The  discharge  rate  through  the  nozzle,  in  this  case  through  the 
throat  region,  can  be  expressed  by  the  continuity  equation 


a)  =  ptvtAt 


(4) 


where 

p  -  density  of  the  gas, 

A  -  cross-sectional  area,  and 
v  -  velocity  at  the  throat. 

The  energy  equation  is 


vt  /  Tt  \ 

-p  =  cp(Tc  -  Tt)  =  cp  Tc(l  -  y~) 


where  the  velocity  in  the  chamber  is  assumed  to  be  negligible. 
Substituting  the  specific  heat  at  constant  pressure 


c  = 


P  V  -  1 


R 


and  introducing  the  isentropic  expansion 

Y-l 


the  velocity  at  the  throat  becomes 


(5) 


(6) 


(?) 


3 


(8) 


The  mass  flow  rate  in  Equation  4  can  be  changed  with 

l 

Pt  /  Pt  \  Y 
Pc  "  ^Pc) 


(9) 


Substitution  of  the  density  of  the  gas  in  the  chamber  with  the  aid  of 
the  perfect  gas  law 


Pc 


(H) 


results  in 


«  =  PcAt 


(12) 


For  convenience  the  right  hand  term  in  Equation  12  will  be  called  a  function 


•-©Vb  [■-(£> 


Yj_ 

Y 


(13) 


which  is  really  the  discharge  function,  depending  only  on  the  ratio  of 
specific  heats  and  the  pressure  ratio  between  chamber  and  throat.  This 
discharge  function  has  a  maximum  value  at  the  critical  pressure  ratio  when 
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sonic  velocity  exists  in  the  throat.  In  order  to  obtain  the  maximum  of  the 
discharge  function,  i.  e.  ,  finding  the  critical  pressure  ratio,  the  derivative 
of  4»  with  respect  to  the  pressure  ratio  has  to  be  taken  and  set  equal  to 
zero. 


2  y+  1  ,  2-Y  J_ 

-(  v  Y  .Ift-E-V  -  (■£-)  Y  Y  -  Hi  (£-Y  1 

d(p/pc)  M  -  1/  2 1_\  pc  /  \pc/  J  L  Y  'Pc'  Y  Vpc/  J 


(14) 


Setting  d'l'/dfp/pc)  =  0  yields  the  critical  pressure  ratio 


jLzX  i 

|  Y  (2/y)  (p/pc)  Y  -  [(V  +  D/V  ]  (p/pc)V 


Tn-; 


(p/Pc)  -  (p/pc) 


(15) 


Equation  15  will  be  zero,  if 


or 


2-Y 

(£-)  y 

Y  + 

Vpc/ 

Y 

^Pc/ 

JL  - 

Y 

-V-1 

Pc 

Vy  + 

1/ 

=  0 


(16) 


(17) 


This  is  the  critical  pressure  ratio  at  which  sonic  velocity  just  exists  in  the 
throat  and  the  discharge  function  has  reached  a  maximum  at 


-  (rhf  Mrr1  •  (ttt) 


max 


Y+_l 

,V  -1 


(18) 


or 


+  =  (-*-Y(  2  y-1 

fmax  \y  +  1/  \V  +  \) 

The  result  shows  that  the  discharge  function  depends  only  on  the  ratio  of 
the  specific  heats  in  the  case  of  critical  conditions. 

For  the  subcritical  case,  the  discharge  function  remains 


1 

— —  ✓ 

Y-l 

+  = 

/Pa\vf  v  f 

(~\  Y 

U 

^pc/  -  iL 

Vpc/ 

J 

With  Equations  3,  12, 

and  13,  Equation 

1  can  be  written 

^=AfP£aPcn-PcAt+(i|-)i  .  (El) 

The  derivative  of  the  gas  mass  in  the  free  volume  of  the  chamber 
can  be  expressed  as 


dp  V 
dt 


(22) 


where  V  is  the  instantaneous  volume  in  the  combustion  chamber  and  p  is 
the  instantaneous  density  of  the  gas.  The  free  volume  is  a  function  of 
time  due  to  combustion  and  consumption  of  the  fuel 

t 

V  =  vc  +jAfaPcndt  (23) 

0 


where 
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m0 

t 

Ua  pcn  dt 
0 


the  VOluiXVft  1TI  t.K p  r*Ka rvvKor.  *»«■  k.. _ 

- - - - UV  WUi  HUU*| 

the  initial  mass  of  the  propellant,  and 


the  volume  of  the  consumed  propellant  or  the 
created  volume  by  combustion  of  propellant 
during  the  time  after  ignition. 


Since  the  chamber  volume  and  the  initial  mass  of  the  propellant 

are  constant,  the  derivative  of  the  instantaneous  volume  with  respect  to 
time  will  be 


dV 

dt 


-  A 


n 


faPc 


Substituting  into  Equation  22  yields 


(24) 


dp  V 
dt 


t 

+  I  AfaPcndt)f  + 
0 


pcAf  ap 


n 

c 


(25) 


The  experiment  shows  that  the  isobaric  combustion  temperature  changes 

only  slightly  with  combustion  pressure,  i.e.  ,  the  temperature  change  is 
negligible. 

Applying  the  equation  of  state  and  taking  the  derivative  of  the  density 
with  respect  to  time  gives 


dPc  I  dPc 

dt  RTC  dt 


With  the  above  relation  the  expression  in  Equation  25  becomes 


dp  V 
dt 


( 


Vc  “ 


m0 

P  f 


J  AfaPcn  dt) 


1 

RT, 


dPc 

dt 


Pc  Af  apc 


n 


0 


(26) 


(27) 
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The  law  of  conservation  of  mass  described  in  Equation  21  can  be 


D 


) 


written  with  the  derived  expressions  as 


( 


Vc  -  —  + 
Pf 


J  Afapcn  dt^ 


1  dPc 


PcAfaPcn  =  AfPfaPcn  ‘  PcAtKRl“)2 


RT„  dt 


0 


(28) 


(Vc  -^  + J  Afapcndt)  =RTcAfPfaPcn(l  -  pcAt*(2RTc)"  .  (29) 

0 


In  this  relation  the  density  in  the  chamber  pc  is  less  than  one  or  two  per¬ 
cent  of  the  propellant  density  p£,  thus  it  can  be  neglected. 


t 

(Vc  -^+jAfapcndt)^-  =  RTcAfp£apcn-PcAttH2RTc)i  .  (30) 

0 

The  above  differential  equation  describes  the  development  of  the  chamber 
pressure  as  a  function  of  time.  The  variables  are 

pc  -  chamber  pressure, 

A^  -  surface  area  of  the  flame  front,  and 

a,  n  -  coefficients  vary  with  operating  pressure  and  pro¬ 

pellant  temperature. 

As  mentioned  previously  the  chamber  temperature  can  be  set  approximately 
constant  as  can  the  gas  constant  or  molecular  weight  of  the  gas  and  the  ratio 
of  the  specific  heats. 

Equation  30  must  be  solved  numerically  on  a  digital  computer.  De¬ 
pending  on  the  design  condition  of  a  solid  propellant  rocket  motor,  different 
cases  have  to  be  distinguished  for  determining  the  chamber  pressure  buildup. 
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If  the  nozzle  is  closed  for  a  moment  during  ignition  by  a  membrane 
or  diaphragm  and  no  efflux  through  the  nozzle  takes  place,  the  following 
relation  has  to  be  applied 


(vc-7F  +  IA£ap"dt)TT  =RTcAfpfaPcn  •  (31) 

0 

This  equation  can  be  used  both  for  the  igniter  propellant  and  for  the  main 
propellant;  only  the  parameters  for  the  propellant  and  the  boundary  con¬ 
ditions  have  to  be  changed. 

At  a  certain  pressure  the  membrane  bursts  and  the  efflux  begins. 
Quasisteady  flow  will  be  assumed,  i.  e.  ,  the  discharge  rate  is  determined 
at  any  given  moment  by  the  equations  for  steady  flow.  Here,  this  equation 
is  valid 


*2 

(Vc  -  ^  +  f  A£aP(>)  =  RTC  AfPfaPc"  -  <2RTC>*  (32) 

0 


Usually,  when  the  diaphragms  burst,  the  chamber  pressure  greatly  exceeds 
the  critical  pressure,  i.  e.  , 


Pa 

Pc 


< 


(33) 


so  the  discharge  function,  which  is  constant,  can  be  introduced. 


Equations  32  and  1  9  yield 


^max 


1 


/  2 

AY-1  /  V  V 

\y  +  1 

!/  \Y  +  1  / 
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Open  End  Nozzle 

The  nozzle  is  open  and  efflux  through  the  nozzle  takes  place  from 
beginning  of  ignition.  Equation  20  gives  the  discharge  function  until  the 
critical  pressure  ratio  is  reached 


•-©V-rHS)’’1 


Then  the  differential  equation  will  be  for  pa/pc  ^  [2/(y+  1)]^  * 


(vc-^^A£apcndt)^=RTcAJPfap 

0 


n 


(35) 


and  for  the  supercritical  phase 


£a 

Pc 


2  \  Y- 1 


ti  Y+  1 

(Vc-^  +  J  AfaPcndt)-^  =  RTcAfPfaPcn  -  (Y)2  (y-|^)2(Y_1)  PcAt(RTc)z  • 


(36) 
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The  time  required  to  build  up  the  chamber  pressure  is  the  sum  of  time  tj 

V 

_  V—  i 

during  which  the  pressure  increases  from  pa  to  pa/[Z/(Y  +  1)1 '  and  the 

V 

Y—  1 

time  t2  during  which  the  chamber  pressure  varies  from  pa/[  2 /(Y  +  1)] 
to  pc.  The  efflux  during  the  first  interval  of  time  is  subcritical  and  ^  depends 
on  pa/p;  the  efflux  during  the  second  interval  is  supercritical  and  4*  is  in¬ 
dependent  of  pa/p. 


Approximate  Solution  of  the  Differential  Equations 

An  integration  of  Equation  30  can  be  accomplished  on  the  basis  of 
assumptions  which  will  follow  later 

*  do  I 

(vc  "  “  +  J  AfaPcn  dt)  ~  RTc  AfPfaPc  ~  Pc  At^  (2RTc)2 

0 

A  new  independent  variable  will  be  introduced  for  the  relative  consumed 
propellant  mass 

t 

4>  =  ~~  J  Afapcn dt  .  (37) 

0 

Equation  30  becomes  with  the  above 

vc  -77  (»  -  =RTcAfpfapcn  -  pcAt4M2RTc)'  •  (38) 

The  derivative  of  the  chamber  pressure  can  be  changed  to 


dpc  dpc  d4>  _  /  dpc  \  /  d<j>\ 
dt  dt  d$  \  d«J>  /  \  dt  / 


and  with  Equation  37 


=  1 L 

dt  mo 


A 


n 

faPc 


(39) 


(40) 
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to 


d?c  Pf  n  QPc 

dt  Sm0Af£pc  d4» 


(41) 


Substitution  into  Equation  38  yields 


vc  -  ~  (1  -  «.) 

c  Pf 


pf 

m0 


Af  a  pc 


n  ape 


n 


d4* 


=  RTC  Af  pfapc*  -  pcAt»J*  {2R Tc) 


{42} 


or  simplified 


( 


VcPf 

m0 


-  1  + 


\  dPc  At  4*  1  -n  _  , 

V  d«|>  “RTcPf~AfaPc  {2RTc* 


Separation  of  the  variable  leads  to 


dp. 


RTcPf-  (At4'/Afa)pc1"n{2RTc)2 


d4> 


[{Vcpf/m0)  -  1  + 


{43) 


The  boundary  conditions  are 


at  t  =0, 


4>  =  0 
Pc  =  Po 


where  pQ  is  the  ignition  pressure,  and 


4*  -  4* 
Pc  =Pc 


The  integration  can  be  carried  out 


Pc 


I 

Po 


dPc 


4* 


1  -n 


RTcpf  -  (At  +  /Ara)  pc  "  (2RTC)2  *  [(Vcpf/m0)  -1+4*1 


rf 


d4> 


(44) 
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The  right  hand  side  integral  can  be  solved  very  easily 


) 


!> 


*> 


4> 


1 

o 


d4>  / 

t(VoPf/m0)  -  1  +  $f  "  V 


VcPf 

m0 


-  1  + 


+) 


♦ 


r(Vc  Pf/m0)  -  1  +  $-j 

1  {VcP£/mo)  -  1  I 


=  in  1  + 


4> 


-  l_ 


(Vc  pf/m0) 

The  expression  for  the  independent  variable  can  be  introduced  again 


I,  =  in 


V 

(p£/mc)  J  A£  a pcn  dt 


1  + 


0 


(Vc  pf/m0 )  -  1 


(45) 


(46) 


or  simplified 


I,  =  in 


1  + 


J  Afapcn  dt 

_0 _ 

vc  -  (m0  /pf)  J 


(47) 


If  it  is  assiuned  that  the  area  of  the  flame  front  is  constant,  which 
is  used  in  many  motors,  and  an  average  burning  rate  is  defined  between 
the  ignition  pressure  pQ  and  the  steady  state  pressure  pc, 


'  m 


/ Po  +  Pc \n 
=  a(~ 2— ) 


(48) 


the  integral  becomes 


Ii 


f  Afa[(Po  +  pc)/2]n 
inV  +  Vc  -  (m0  /Pf)  / 


(49) 
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For  the  solution  of  the  left  hand  side  integral  in  Equation  44  some 
assumptions  nave  to  be  made  such  as  constant  area  of  flame  front  and 
constant  discharge  function  in  the  case  of  supercritical  condition,  which 
is  true  when  the  ignition  pressure  pQ  is  above  the  critical  pressure.  The 
combustion  index  of  many  propellant  types  are  in  the  neighborhood  of 
n  =  2/3,  so  it  will  be  used  for  the  integration. 


Pc 


I 


dp. 


Pc 


■I 


dp. 


=  I? 


p  RTcpf  -  (At4*/Afa)  pc1-n  (2RTc)1/2  p0  C-dpc  h 


(50) 


where 


C  =  RTcPl 


and 


At^ 

d  =  A-r  <2RTc» 


which,  by  dividing  out  the  integral,  becomes 


Pc 


'2  /  j2  »  -^3 


12  =I  + 


C  (C  /d  )pc 


C  d2  C 


C-dpc 


Vs  j 


|dp. 


Integration  of  the  terms  results  in 


t  -  3  2 L 

lz  "2d  Pc 


3C  i/3 
"d1  Pc 


£l»,C-dpc*>] 


Po 


and  setting  the  limits  yields 


(51) 


(52) 
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2  C  -  dp  1/3 


T  2_,  2/3  ft,.  3C  .  1/3  1/3.  3C2  _ 

Iz  “  2d  *P°  ”  Pc  *  +  "cF  (Po  '  Pc  *  *  IF* in  c  dp  i/?> 


3C2  T  1  .  ft  1/3 x  .  d  .  2/3  ft.  1  , 

d2  [C  "  Pc  *  +  2C2  Po  “  Pc  ~  d  in 


C  -  dp 


.  (53) 


C  -d 


The  whole  expression  of  the  solved  integral  I2  is  given  with  replacement 
of  the  constants 


*P/**V*TC/  !  ,  *  *  ,  At+'2RTc)'/2  * 

24'zAt2  VPfRTc  2A£apf2R2Tc2 


'  Pc.  ) 


Af  a 


Pf RT  -  [At  +  (2RTc)‘/2/Afa]pc‘A 


1/Z  I  K  1  ^3 


A>{2R.Tcr  UfRTc  -  [At«M2RTcr /Afa]Po 


The  solution  of  Equation  44  permits  the  determination  of  the 
chamber  pressure  as  a  function  of  time  which  is 


3pf2  a2  Af2  RTC 


ZV  At‘ 


1  ,  ./3  4„  At"'2RTc)'2  ,  *  * 

— —  (Po  -  Pc  )  +  -  VPo  -  Pc  > 

2Afapf2  R2  Tc2 


)*■  A.1  VPfRTC  ?Ara  nsZ  R2  T  2 


At^(2RTc) 


Pf RTC  -  [At4>(2RTc)1/z/Afa]  pcV3' 
Pf  RTC  -  [  At^(2RTc)l/?  / Af  a]  pQl/3. 


r  Af  a[(p0  +  Pc)/2]  t 

in  l  +  Vc  -  (m0  /pfj 
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If  the  area  of  the  flame  front  is  a  variable,  i.  e. ,  a 
and  the  combustion  index  different  from  n  s  2/3  especially 
sure  ranges.  Equation  30  has  to  be  integrated  numerically 
an  accurate  chamber  pressure  curve. 


function  of  time 
in  certain  pres- 
for  obtaining 


3) 


4) 


3) 
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CHAMBER  PRESSURE  DECAY 


Decay  During  Supercritical  Phase 

After  burnout  the  combustion  chamber  can  be  considered  as  a  gas- 
filJ.ed  vessel,  in  which  pc,  pc  and  Tc  exist  as  initial  conditions.  It  is  as¬ 
sumed  that  the  discharge  of  the  gas  through  the  nozzle  is  approximately 
quasistationary,  i.  e.,  the  discharge  rate  is  determined  at  any  given 
moment  by  the  equations  for  steady  flow. 

The  law  of  conservation  of  mass  applied  to  the  gasecus  mass 
contained  in  the  combustion  chamber  can  b<'  introduced  as 


dp  V 
dt 


u> 


(56) 


where  «  is  the  discharge  flow  rate  out  of  the  chamber. 

Since  the  volume  of  the  chamber  does  not  vary  with  time,  Equation 
56  becomes 


v  dp.  * 

vc  dt  "  “ 


(57) 


Using  the  derived  relation  for  the  discharge  rate  in  Equation  10 
the  expression  will  be 


dp  _ 


£  /pt\v 


c  dt 


-  pAt<2RT}2  (~) 


Ill 

y  r,  _  v 

Y  -  1  L  v  P  / 


] 


(58) 


where  the  dis  '.harge  function  ^  is  introduced  again  by 


Y"  1 
/Pt\  y 


-C'.  i 


(59) 
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I 


Dividing  by  the  initial  chamber  condition  yields 


±Z.  =.£- 


dt  p 


1  i 

Aj^RTc)1  (^)2  4.  . 

c 


'C  Pc 

Following  the  adiabatic  relation  and  arranging  gives 


1 


A 

dt 


(JL^ 

^ } 
II 

U 

— «v 

Wc> 

'  Vc\ 

Y+l 


The  derivative  can  be  written  as 

1 


i/i-V  =i('-£.y'  'ifjE-N 

it\pc/  Y\p  c)  dtVpr/ 


1  -Y 


=iYJ£.\ v  i/p.'j 

Y\pc/  dt\pc/ 
Substituting  into  Equation  61  yields 


3Y-1 

=  .^ly  (JL)2y 

dtVp  J  Vc 


(f~)2Y  ^(2RTc) 

X  Pc-' 


(60) 


(61) 


(62) 


(63) 


This  first  order  differential  equation  can  be  solved  by  separation 
of  the  variables  as  follows 


d(p/pc) 


3Y-I 

<Mp/pc)  Zy 

The  boundary  conditions  are 


=  ___Y(2RTc)2  dt 


(64) 


t  -  0  -  =  1 

Pc 


t  =  tj 


_P_  _  JL 
Pc  "  Pc 
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F.nnatmn  Ai 1  ir>t  p/t  rooulfo  i  r' 

- x - - - -fe.— *  ' —  '**vw  1A* 


p/pc 

r  ,d(p/fc)  = 

J  3V-1 

i  +<p/Pcr*r 


Y(2RTc)2  t, 


Since  the  discharge  function  ^  depends  on  the  pressure  ratio,  two 
cases  have  to  be  distinguished,  the  supercritical  and  the  subcritical  con¬ 
ditions.  For  the  supercritical  condition  the  pressure  ratio  is  as  developed 


previously 


V 

2  YY-1 


<  (_i_Y 

P  -  Vy  +  i/ 


and  the  discharge  function  has  reached  a  maximum  with  sonic  velocity  in 
the  throat 


+  =  (^-f  (—  -V~‘ 

max  \Y  +  1/  \V  +  1  / 


Equation  65  becomes 


At  —  C 

—  Y(2RTc)2t1  =  -  j 


P/Pc 


(P/Pr) 


3Y_1 

2Y 


_d(^L)  . 

1  VPC/ 


1  [Y/CY  +  1)]2[2/(V  +  1)]V“! 


The  integral  will  be  calculated  separately 


p/pc  31M 

I  (£)  ”  *(fc)  ■ 


3Y-1 


- 1 - W 

(3Y  -  1)/2Y  +  1  Pc/ 


(£) 
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(67) 


P/Pc  Itl 
f  (JL\  2Y 
J  VPc^ 


Y-l 

d/_E_\  =  P_\“  2V 

\pc/  V  -  1  \pc/ 


Y-l 


With  this  result  Equation  66  becomes 


At  — 

—  V(2RTc)2t, 

vc 


_ 1 _ 

1 

['v  /( Y  +  1)]^[2/(Y  +  l)]7"1 


•  (68) 


Rearranging  the  above  equation  and  expressing  the  chamber  pres¬ 
sure  explicitly  as  a  function  of  time  gives 


_P_ 

Pc 


At  A 

1  +  —  Y^RT^t, 
vc 


1 

Y  -  1  /  Y  \I  /  2  \  V-l 

2Y  VY  +  1/  \Y  +  1/ 


_2Y 

Y-l 


or  simplified 


_P_ 

Pc 


Y+l 

Y-l/  2  \2(Y-1) 

2  \Y  +  l) 


(YRTC)2 


2Y 

Y-l 


(69) 


Equation  69  shows  that  the  efflux  at  supercritical  condition  is  independent 
of  pa/pc. 


Decay  During  Subcritical  Phase 


For  this  phase,  Equation  65  is  also  valid  which  is 


Av  1 

~  Y(2RTc)2t 


3Y-1 


1 


(70) 
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With  the  subcritical  discharge  function 


L'-(r) 


Pa\  V  U 


Taking  into  account  the  initial  condition 


t  =  0 


—  =  ( 

p  \V  +  1  / 


2 _ \  y-1  _  Pa  Pc 

+  l)  "P  Pc 


-£-  = 
Pc 


pj2/(V  +  D] 


Pc  Pc 


the  relation  becomes 


~y(2RTc)2t2  =  - 
v  c 


p/f 

/ 


pa/  pc[2/(Y+l)] 


_  3Y~1 

(p/p  )  2Y  d{p/p  ) 


[Y/(Y  -  l)]2  (Pa/P)Y  [1  "  (Pa/p)  ^  3 


For  better  integrability  some  minor  changes  must  be  performed 
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ni>- 


YAtt2{[2Y/(Y  -  1)]RTC}2 
— - 

_  3^~l> 

(p/pc)  ^  d(p/Pc) 

3EE  V-l 

U  “  (Pa/Pc)  Y  (p/Pc)  Y  ]2 

(74) 

The  above  integral  can  be  calculated  only  by  the  numerical  method. 
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THRUST  BUILDuV  ANjd  ujlCAY 

Thrust  Buildup  and  Decay  During  Supercritical  Phase 

With  the  aid  of  the  previously  determined  derivations  the  chamber 
pressure  as  a  function  of  time  can  be  calculated  for  the  buildup  and  the 
decay,  respectively.  Applying  the  law  of  momentum  on  a  rocket  nozzle, 
the  thrust  is  given  by 


where  u>  is  the  mass  flow  rate,  vg  is  the  exit  velocity  or  the  end  velocity 
where  separation  of  the  gas  from  the  wall  occurs,  Ae  is  the  area  in  the 
nozzle  at  separation  and  pe  is  the  static  pressure  at  this  point. 

In  the  above  equation  for  the  thrust  the  correction  factor  X.  will  be 
assumed  to  be  unity. 


1  t  cos 

X  = - - - *  1  (76) 

where  is  the  half  nozzle  angle  for  the  divergent  part  of  a  conical  nozzle. 

The  mass  flow  rate  through  the  nozzle  is  given  by  Equation  12 
which  is 


23 


or 


w  =  pc  At 


u> 


4> 


(77) 


where  ^  is  the  discharge  function. 

For  the  critical  phase  where 


Pa 

Pc 


< 


the  discharge  function  reaches  a  maximum  and  is  a  constant 


4» 

max 


In  this  case  the  mass  flow  becomes 


1 


Cu 


=  Pc  Why  (yhr 


(78) 


or 


10  = 


-  PcAt(RT  )2  (  y  +  1  ) 


Y+  1 
2  \2(V-1) 


(79) 


From  the  energy  equation 


Cp(Tc 


(80) 


and  the  equation  for  isentropic  expansion 
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(81) 


Y-i 


the  end  velocity  ve  follows 


For  obtaining  the  area  ratio  where  the  separation  occurs  the  con¬ 
tinuity  equation  must  be  applied. 


Pe  ve  Ae  -  PtvtAt 


(83) 


With  the  end  velocity  found  in  Equation  82  and  the  velocity  of  sound  in  the 
throat,  the  expression  becomes 


Y-l 


2Y 


fl  -  Y 


RTJ1  - 


lY  -  1  CL  \Pc/ 


2  Pt  ?■ 

=  pI  (vRTt)  At 


"tit  (VRTc)20f)2  At 

Pc  Pe  ' *c/ 

1 

1 

_/  2  V"1 

~  \Y  +  1/ 

l  l 

=  (Y)2  (RTC)2 

/Pc\V  /  2  \2 

\pe/  VY  +  1/ 

the  area  ratio  is  found  by 

Y+  1 

(Y)2  [2/(Y  +  l)j2(Y'1) 

1 

Y-l 
'  1 

(84) 


Y+  1 


(Pe/Pc)  {[2Y/(Y  -  l)j[l  -  (pe/pc)  ]}: 


(85) 
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The  thrust  is  found  by  substituting  into  Equation  (5: 


F  =pcAt(Y)^  (y^-j 


Y+  1 

Ill 

\2{Y-1) 

{  2V 

•  O '  ]} 

) 

\Y-> 

At(Y)H2/(Y  +  1)] 


v+i 

2(V-1) 


(86) 


Y-l 


(Pe  ‘  Pa) 


(Pe/Pc)  {[2Y/(Y  -  1)][1  -  (Pe/Pc)  ]) 


This  equation  is  valid  for 


Y 


At  the  design  point  of  the  nozzle  the  gas  is  expanded  down  to  the  nozzle  end 
where  the  exit  pressure  is  equal  to  the  ambient  pressure  pa.  If  the  chamber 
pressure  is  lowered,  the  gas  is  overexpanded  and  the  exit  pressure  can  be 
much  lower  than  pa;  even  separation  from  the  wall  occurs  and  the  separa¬ 
tion  point  travels  upstream  with  decreasing  chamber  pressure.  During 
overexpansion  oblique  shock  waves  exist  in  the  nozzle. 

From  the  standpoint  of  thrust,  the  nozzle  may  well  be  cut  off  at 
the  separation  station,  since  the  internal  and  external  pressures  are  nearly 
in  balance  beyond  this  point.  The  thrust  of  a  rocket  motor  with  separation 
in  the  nozzle  is  calculated  on  the  assumption  that  the  nozzle  area  ratio  is 
not  that  corresponding  to  the  actual  exit  but  rather  that  of  the  separation 
station.  Therefore,  it  is  necessary  to  be  able  to  predict  the  location 
of  separation. 

Experimental  results  on  rocket  motors  which  vary  considerably 
from  author  to  author  will  be  used  for  finding  the  exit  pressure  where 
separation  of  the  gas  from  the  nozzle  wall  occurs.  Summerfield  recom¬ 
mends  an  average  exit  pressure  ratio  of  pe/pa  =0,4  throughout  all  the 
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expansion  ratios,  but  this  value  i»  not  realistic,  especially  when  the 
separation  exists  in  the  neighborhood  of  the  throat  or  when  only  sonic 
velocity  is  the  maximum  velocity  due  to  the  available  pressure 
ratio.  In  Figure  1  the  exit  pressure  can  be  found  as  a  function  of  the 
chamber  pressure  ratio.  With  this  exit  pressure  pe  the  thrust  in  Equation 
86  is  easily  calculated. 

Thrust  Buildup  and  Decay  During  Subcritical  Phase 
In  the  case  of  subcritical  pressure  ratio 


Pc  -  \Y  +  1/ 


2  \  Y-l 


the  pressure  at  the  throat  is  equal  to  the  ambient  pressure 


Pe  =Pa 


Using  Equation  77  the  mass  flow  rate  becomes 


(87) 


and  with  Equation  82  the  end  velocity  is 


Y-l 


* 


trtc 


1  - 


/  Pa\  Y  « 

v*)  J 


1 

2 


(88) 


Finally,  the  thrust  cam  be  expressed  as  follows  for  the  subcritical  phase 


[■  ■  & 


(89) 
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Figure  1.  Exit  Pressure  at  Separation  as  a  Function  of  Chamber  Pressure  Ratio 


RESULTS 


Expressions  for  the  chamber  pressure  and  thrust  have  been  derived 
as  a  function  of  time  for  the  supercritical  and  subcritical  conditions  during 
buildup  and  decay.  Special  attention  was  given  to  certain  rocket  motor 
design  conditions,  such  as  open  end  nozzles  and  closed  nozzles  by  diaphragms 
during  ignition.  The  results  are  differential  equations  which  must  be  solved 
numerically  for  accurate  purposes  and  in  the  case  of  variable  parameters 
in  the  equations.  In  addition  closed  form  solutions  have  been  found  based  on 
assumptions  such  as  constant  propellant  burning  area  and  average  burning 


CONCLUSIONS 


In  miny  cases  the  total  impulse  during  the  thrust  buildup  and 
decay  period  is  neglected  by  calculating  average  performance  values. 

In  order  to  evaluate  the  exact  total  impulse  the  transient  conditions  have 
to  be  taken  into  account.  The  instantaneous  chamber  pressure  and  thrust 
which  are  affected  by  the  characteristics  of  the  propellant  and  the  design 
of  the  motor  must  be  determined  throughout  the  operating  time  of  the  rocket 
motor. 
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